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Abstract 



We construct isospectral partner potentials of a complex VT- invariant potential, viz., V(x) = 
—Vi sech 2 x — %Vi seek x tank x using Darboux's method. One set of isospectral potentials 
are obtained which can be termed 'Satellite potentials', in the sense that they are of the same 
form as the original potential. In a particular case, the supersymmetric partner potential 
has the same spectrum, including the zero energy ground state, a fact which cannot occur in 
conventional supersymmetric quantum mechanics with real potential. An explicit example of 
a non-trivial set of isospectral potential is also obtained. 



Key words : VT- symmetry, non-Hermitian Hamiltonian, isospectral potential, Darboux's 
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I. Introduction 



Quantum systems characterized by non-Hermitian Hamiltonians have been studied from time 
to time, because of their applications in scattering problems [1]. For this reason, complex 
potentials were also termed as optical or average nuclear potentials. The interest in this field 
has been revived fairly recently after it was conjectured by Besis, Bender and Boettcher, and 
others, that Hermiticity of the Hamiltonian is not essential for a real spectrum [2,3]. Several 
non-Hermitian PT-invariant complex potentials have been found to possess a real spectrum, 
though the corresponding Hamiltonians are non-Hermitian [4,5]. Moreover, PT-invariant mod- 
els admit some of the properties of the usual Hermitian ones, viz., super symmetry, potential 
algebra, quasi-solvability, etc. [6-17]. Such non-Hermitian PT-invariant Hamiltonians have 
found applications in many areas of theoretical physics — nuclear physics, field theories when 
studying Lee- Yang zeros, localization-delocalization transitions in superconductors, theoretical 
description of defraction of atoms by standing light waves, and also the study of solitons on a 
complex Toda lattice [7]. 

In this study, our aim is to construct new complex potentials via the Darboux's method 
[7-10], for which the corresponding eigenvalue problem can be solved exactly. The constructed 
potentials are not necessarily VT- invariant, but still give rise to a real and discrete spectrum, 
provided the original potential admits real energies only. 

The organization of the present note is as follows. To make it self-contained, we give a brief 
review of PT-symmetry in Section II. In Section III, we give an outline of the Darboux's method 
for constructing isospectral partner potentials of a given known PT-symmetric potential. In 
Section IV, we illustrate our approach with the help of an explicit example, viz., the one- 
dimensional PT-invariant potential 

V(x) = —Vi sech 2 x — iV 2 sech x tank x V\ > (1) 

Section V is kept for conclusions and discussions. 



II. VT- symmetry 

The Hamiltonian H for a particle of mass m, in a complex potential 
V(x) = V R (x) + iVj(x) is given by 

H is said to be PT-symmetric when 

VTH = HVT (3) 
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Here V is the Parity operator acting as spatial reflection, and T stands for Time Reversal, acting 
as the complex conjugation operator. Their action on the position and momentum operators 
are given by : 

V : x — > — x, p — > — p 
T : x — > p — > — p, i — > — i 
Hence, in explicit form, the condition for a potential to be PT-symmetric is 

V*(-a;) = V{x) (4) 

The commutation relation 

[x,p]=ih (5) 

remains invariant under VT for both real as well as complex x and p. 

It is worth mentioning here that though the Hermiticity of the Hamiltonian may be replaced 
by the weaker condition of PT-symmetry, the latter is not sufficient for the reality of the 
spectrum. Various authors have studied several one-dimensional non-Hermitian PT-symmetric 
models and shown that such Hamiltonians exhibit 2 types of behaviour — 

(i) In the unbroken PT-symmetry phase, the eigenfunctions of H are also eigenfunctions of 
VT, and the energy spectrum is real and discrete. 

(ii) In spontaneous breakdown of PT-symmetry, though the potential retains PT-symmetry, 
the corresponding wavefunctions do not, and the energy eigenvalues exist as complex conjugate 
pairs. 



III. Darboux's method 

The Darboux's method [7-10] relates the spectral properties of a pair of standard Schrodinger 
Hamiltonians 

H - = ~^ + v ^ < 6 > 

We assume that the spectral properties of one of these Hamiltonians, say H + , is exactly known. 
Thus 

H + 4> n {x) = E n <f> n (x) (7) 

where the eigenvalues E n and the corresponding eigenfunctions (j> n (x) are known explicitly. The 

spectrum is assumed to be discrete such that E < E\ < E 2 < ■ ■ ■. 

Let there exist a linear operator A, such that it obeys an intertwining relationship 

AH + = H_A (8) 

Hence the functions ip n = A(p n are eigenfunctions of if_ with the same eigenvalues E n 

H_^ n (x) = E n il> n {x) (9) 
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A general form for an intertwining operator A obeying (3) may be given by [7,8] 

A = E/*(*)|* (io) 

where fk(k = 0, 1,2, N — l)are at least twice differential functions and fx is an arbitrary 

constant. 

For simplicity of calculations, we work in units H = 2m = 1. 
We take the simplest non-trivial choice for A(N = 1) 

A=-± + m (id 

Putting (6) in (3) yields 

[V _ _ V+ + 2f] ± - [(V. -V + )f + Vi-f"]A = (12) 

where prime denotes differentiation with respect to x. This relationship (12) puts certain 
restrictions on the functions V+(x), V-(x), and f(x), viz., 

V4x)=V+(x)-2f'(x) (13) 

[V_(x) - V + (x)} f(x) + V[{x) - / "(x) = (14) 
Substituting (14) in (13), and integrating, we obtain 

f\x) + f\x)-V + (x) = -e (15) 

where e is an arbitrary integration constant, in general complex. 
Putting 

- U -M (16) 

(15) can be cast in the form of Schrodinger-like equation 

d 2 u(x) 



dx 2 



+ V^ar) = eu(x) (17) 



e is sometimes called the factorization energy. It is worth noting here that u(x) need not be 
square- integrable. So we are not restricted to normalizable solutions of (17). However, for A to 
be well defined, u must not have any zeroes on the real line. For this, it requires the condition 
e < E . 

In terms of the function f(x), the two potentials are expressed as 

V ± (x) = f 2 (x)±f'(x) + e (18) 
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At this point it is obvious that for e = 0, the two potentials V±(x) are the supersymmetric 
(SUSY) partner potentials, and H± are the SUSY partner Hamiltonians. From (18), the partner 
Hamiltonians may be expressed in terms of the linear operator A as, 

H + = A + A + e (19) 

H = AA + + e (20) 

We shall apply the above approach to construct a new complex potential V-(x), which is 
isospectral with a known PT-symmetric potential V+(x) with real eigenvalues, obeying the 
relationship 

^ 2 Kf} 2 - y+w+2e < 21 > 

Though PT-symmetry may be either restored or broken in the new potential, the partner 
Hamiltonian if_ has the same eigenvalues as H + (with the possible exception of the ground 
state). The corresponding eigenfunctions are given by 

, / n f dcbJx) uf(x) , . .] 

M*) = C n --^ + -l^f <t>n (X) 22 

[ ax u(x) J 

Here c n is the normalization constant, given by 

|c n |~ 2 = (A<p n \A(j) n ) ^ 
— E n — e 

We shall make our ideas clear with the help of an explicit example in the next section. 



IV. Explicit example 

As an explicit example, we shall consider the one-dimensional PT-invariant potential [10,18] 

V(x) = —Vi sech 2 x — 1V2 sech x tank x V\ > (24) 

This non-Hermitian Hamiltonian has the interesting property of having both real and complex 
discrete spectrum, depending on the relative magnitudes of its real and imaginary parts, viz., 
Vi and V 2 . 

(i) Foi\V 2 \<V 1 + l 

"PT-symmetry is unbroken and energies are real. 

(ii) For \V 2 \ >V 1 + \ 

PT-symmetry is spontaneously broken and energies are complex conjugate pairs. 

We shall restrict our attention to the unbroken PT-symmetric phase, with real energies only. 
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To solve the Schrodinger-like equation 

d 2 u(x) 
dx 2 



+ V(x)u(x) = eu(x) 



with V{x) given by (24), we make the following substitutions : 

1 — i sink x 



u (x) = z- p (1 - z)- q X (z) 
Then the differential equation satisfied by x( z ) is 



z(l-z)x " + 



provided 



-2p + - + z (2p + 2q-l) 



1 + s 

v — — ± — 
y 4 2 

9 = -4 ± 2 



X - 



(p + q) 2 -\ 2 X = 



t = J Vl -v 2 + l 



'^1 + ^ + - 



A 2 = -e 

The most general solution of (27) is 

X(z) — a F(a, b, c, z) + pz 1 ^ (f - z) c ' a ' b F (1 - a, 1 - b, 2 - c, z) 

where 

a = —p — q + A 
b = —p — q — A 
1 

c = -2p + - 



(25) 
(26) 

(27) 

(28) 
(29) 

(30) 

(31) 
(32) 

(33) 

(34) 
(35) 

(36) 



So the general form of u is 

u = a z' p (1 - z)" 1 F (a, b, c, z) + (3z p+1/2 (1 - z) q+1/2 F (1 - a, 1 - b, 2 - c, z) (37) 

Since u must not have any real zero, so a must be non-zero. This allows us to put a — 1. 
However, the proper choice of the sign in the expressions for p and q is extremely important. 
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For the potential given in (24), the real energy eigenvalues and the corresponding eigenfunc- 
tions are explicitly given by [18], 



E n = - jn + l - - X - (s + i)} n = 0, 1, 2, .... < 



s + t-1 



where P£' (z) are the Jacobi polynomials expressed as [19] 

P a /{Z) = 2-» ± ( n + ° ) ( ^ ^ ) (* - 1) — (Z + iy 



(38) 



a / \ I n \ f 1 ~ i sinh x \ ( l ±i sinh x \ 2v -\-i q -\ ( . . , . 
Mx) = \ n _2 p -1 )[ 2 J \ 2 J 2 {isinhx) (39) 



(40) 



Case (i) 

First let us consider the simplest case given by 

(3 = 



Hence 



Using the relationship [20] 
u(x) reduces to the simple form 

u{x) = a 



b = c 

\ = p-q-}- 
e = -X 2 = -[p-q-l/2} 2 

F(a,b,b,z) = (l-z)- a 
( 1 — i sinh x\ p f 1 + i sinh x 



9+1/2 



Thus f(x) = u'(x)/u(x) is obtained to be 



f(x) = i 



p+ [q+- 



sech x — 



P-\Q+2 



tanh x 



The isospectral partner potential of (24) as obtained from (21) turns out to be 



V-(x) 



v 1 



sech x 



(41) 
(42) 

(43) 
(44) 



(45) 



— I 



sech x tank x 



(46) 



with the corresponding eigenfunction (from (22)) ip n as 

, 1 — i sink x\ p ( 1 + i sink x\ 9 



i J (iq + [i sech x + tank x] P n 2p 2 ' 2<? 2 (i sink x) 



n-2p-± 



(p + q) 2 -E n 



! ] (iq + -\ [i sech x + tank x] P n 2 ± + 2 ' 2q+ 2 (i sink x) 



-2p + l/2 \n-2p- 

(47) 

This is the special case of the so-called satellite potentials, with V+(x) and V-(x) having 
exactly the same form, but with different coefficients. 
We illustrate this with the help of some simple examples. 
Let V 1 = 25, V 2 = 5 
Thus 

V+(x) = —25 sech 2 x — hi sech x tank x (48) 

with energies 

£ n = -{n + I- 5 } , n = 0,l,2,3,4 (49) 

Then using equations (28) - (31), the values of p and q are 
p = \ or -3 

q = 2 or -| 

In order that the wave function vanishes asymptotically, only the positive sign of the discrim- 
inant is allowed for p. Therefore, p = 5/2. The two values of q give rise to two interesting 
cases. 



a) q = -5/2 

f(x) for this particular case becomes 



11 i 
fix) = — tank x — -sech x 
J v ' 2 2 



and the isospectral partner potential of (48) turns out to be 

V-(x) = —16 sech 2 x — Ai sech x tanh x 

with energies 

£ n = -{n+I-4} , n = 0,l,2,3, 

It is easy to check that V±(x) share exactly the same energy spectrum with the exception of 
the ground state, as shown in Table 1 below. 



(50) 

(51) 
(52) 
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Table 1 



n 


E n [V + (x)} 


E n [V-(x)] 





-81/4 


-49/4 


1 


-49/4 


-25/4 


2 


-25/4 


-9/4 


3 


-9/4 


-1/4 


4 


-1/4 





b) q = 2 

This is the special case when p = q + 1/2 , so that e = 
and f(x) becomes purely imaginary. 

f(x) = i 5 sech x 

The isospectral partner reduces to 

V-(x) = —25 sech 2 x + 5 i sech x tanh x 

Thus 



(53) 



(54) 



V ± {x) = f{x)±f{x) (55) 

The partner potentials (48) and (54) are totally degenerate. They share identical energies, 
including the zero energy ground state, and supersymmetry (SUSY) is broken. This scenario is 
quite different from conventional SUSY breaking. In the conventional supersymmetric quantum 
mechanics (SUSYQM), SUSY is broken when the zero state energy does not exist and both 
V±(x) have the same spectrum. The ground state wavefunctions <f>o(x) and ipo(x) of the partner 
potentials (48) and (54), are respectively given by 



M x ) = 



i sinh x 



1 + i sinh x^ 



sech x 



sech xyl + cosh x + 



tanh 



a/1 + cosh x 



(56) 



ip (x) = 



9 
2 

sech 5 x 



dx ' u(x 
[tanh x + i sech x] (f>o(x) 

tanh x VT+ cosh x 



(57) 



+ i sech x y 
apart from normalization constants. 



sech 2 x 



VT + cosh x J 



a/1 + cosh x 



+ Vl + cosh x — 



1 



a/1 + cosh x 
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Case (ii) 

Next we consider the particular case 



(3 = 0, b = a+-, c = 2a 
This gives the following values for 

x 1 1 1 

A = -4' £ = -16' q = ~2 

Hence, for this particular case, V 1 = V 2 = V (say) in (24). Using the relationship [20] 
F[a,a + 1 2a, = 2 2ffl - 1 (l - z)- 1 ' 2 [l + (1 - z) 1 ' 2 ] 1_2 ° 



u(x) reduces to 



u(x) = a 2~ 2p - 1 ' 2 



1 — % sink 



1 + i sink x 



1/2' 



2p+l/2 



Thus f(x) = u'(x)/u(x) is obtained to be 



fix) = -tank x — — sech x + i ( 2p H — ^ sech x 
4 4 V 2) 



1 + i sink x 



1/2 



The isospectral partner potential of (24) viz., 

V+(x) = —Vq sech 2 x — iVo sech x tanh x 
is obtained from (21), and turns out to be 



V-(x) = 



1 / 1 
4 + { 2p+ 2 



V n 



sech x — i 



1 / 1 
4 + { 2p+ 2 



sech x tanh 



+ (2p 



sech x + i sech x tanh x 



1 + i sinh x 



1/2 



or more explicitly, 



V-(x) = — Vr(x) sech 2 x — iVj(x) sech x tanh x 



with 



V R ( X ) =a+(p+\) Vco 2 x+1 



cosh 2 x 



V cosh x + 1 

VAx) = a -(p + \)\ V cosh x + 1 H — , 1 
V 4/ [ ^cosh x + 1 J 
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where 



and 



a = V + 2p + - - V 



(66) 
(67) 



It is worth noting here that the coefficients of (sech 2 x) and (sech x tank x) are pure numbers 
in the expression for V+(x), whereas they are functions of x in V-(x). From equation (38) it is 
essential that s and hence V be large for sufficient number of energy levels, as t — 1/4. This 
prevents p from taking the value —1/4. So the isospectral partners V±(x) in equations (61) and 
(63), cannot be reduced to satellite potentials for any physical value of p in this particular case. 
The corresponding eigenfunction is calculated to be 



i/>n(x) 



-> I' ' ~\) S ' X + (/' "J ) ttlllll .)■ T / 



2p + — ) sech x 



1 + % sink x 



1/2- 



n \{l — isinhx\ /l + ism/ix\ 2 n -2 P -ii,. . , , 

c 1 \ r n 2 2 {i sink x) 

n - 2p - 1/2 J \ 2 J \ 2 J y ' 



(p — 1/2) 2 — E n (\ — i sink x\ p (\ + i sink x" 



-2p + 1/2 



n- 1 
n-2p- 1/2 



P n 2p+2 ' 2 (i sink x) 



(68) 



V. Conclusions and Discussions 

To conclude, it is shown in this paper how to generate equivalent SUSY partners of the complex 
VT- invariant potential, viz., 

V[x) = —V\sech 2 x — i V^sech x tank x 

with the help of Darboux's method. The form of the newly constructed potentials depend on 
the choice of f3 and e, but they share the same energy spectrum as the origial potential, with the 
possible exception of the ground state. If the original potential is so chosen that it admits real 
eigenvalues only, then we can obtain a series of non-trivial complex potentials generating the 
same real-valued spectrum. The new wavefunctions are also easily obatined by this approach. 

Depending on particular values of V\ and V2, the isospectral partners may be of similar nature 
with just the coupling constants taking different values. We term these as satellite potentials. 
In the special case p = 5/2, q = 2, e = — (p — q — 1/2) 2 — 0, the SUSY partner potentials have 
the same spectrum including the zero energy ground state, a situation which cannot occur in 
conventional SUSYQM with real potential. 



12 



We have also constructed a non-trivial partner potential, and compared its real and imag- 
inary parts graphically with those of the original potential, in Figures 1 and 2 respectively. 
It is interesting to note that even in the non-trivial case, the partner potentials have similar 
geometric form. 

Further constructions with non-zero (5 and other values of the integration constant e are left 
as a future exercise. 
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Fig. 1 

solid line : Imaginary part of V (x) [eq.(61 )] 
broken line : Imaginary part of V + (x) [eg. (63)] 



-x x 

Fig. 2 

solid line : Real part of V + (x) [eq.(61 )] 
broken line : Real part of V (x) [eq.(63)] 



